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. . . We present a simple method, combining the density-matrix renormahzation-group (DMRG) 

\ algorithm with finite-size scaling, which permits the study of critical behavior in quantum spin 

. chains. Spin moments and dimerization are induced by boundary conditions at the chain ends and 

' these exhibit power-law decay at critical points. Results are presented for the spin-1/2 Heisenberg 

antiferromagnet; an analytic calculation shows that logarithmic corrections to scaling can sometimes 
51^ ■ be avoided. We also examine the spin-1 chain at the critical point separating the Haldane gap 

^ ' and dimerized phases. Exponents for the dimer-dimer and the spin-spin correlation functions are 

consistent with results obtained from bosonization. 
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I. INTRODUCTION 



, Quantum critical points are characterized by fluctuations over all length and time scales and by the appearance of 
power law scaling. In this paper we present a simple but powerful numerical method to access quantum critical points 
in one-dimensional systems. The method combines the density-matrix renormalization-group (DMRG) algorithm and 
■ O finite-size scaling ideas. We illustrate the method by applying it to several well-understood quantum spin chains. In 
C a second paper to follow we, apply the method to new classes of supersymmctric spin chains which describe various 
^ ,: disordered electron systemstJ. P 
I i' The development of the density- matrix renormalization-group (DMRG) algorithm by WhiteH represented an im- 
portant improvement over previous numerical methods for the study of low dimensional lattice models. It has been 
applied to a wide variety of systemstl. The DMRG approach was first used to study the ground state properties 
and low-energy excitations of one-dimensionaLchains. It has been extensively applied to the study of various spin 
chains. Low-lying excited states of the spin-lQ □ and spin-l/2Q Heisenberg antiferromagnets have beei|i|-C|alculated. 
^ Likewise, spin-1 chains with quadratic and biquadratic interactionajtl, a spin-2 antiferrtjmiagpetic chairil30, spin-1/2 
^— I . and spin-1 chainSpWith dimerization and/or frustration (nextrngarpit- neighbor coupling)ll5il3, and frustrated spin-3/2 
' and spin-2 chainalJ have all been studied. Edge excitationsEflljIljJia at the ends of finite spin chains and the effects of 
, perturbations such as a weak magnetic field coupled to a fnsa. sitesc3 have been considered. j-B,andomness in the form of 
' random transverse magnetic, field in a spin-1/2 XY modelEil, random exchange couplipffSLJ, and random modulation 
patterns of t he-exchanger ^Ir ^, has been examined. Finally, alternating spin magnitudescJ, the presence of a constantcS 
or a staggeredEZi magnetic-field in a spin-1 chain, bond dopin^a, the effects of a local impuritycj, and interactions 
' with quantum phononsE3til have also been considered. 
• Most of the above work involves systems in which the first excited state is separated from the ground state by 
' a non-zero energy gap as the DMRG works best for gapped systems. First attempts to extract critical behavior 
Q , of gapless systems used the DMRG to generate renormalization transformations of the coupling constants in the 
O ■ Hamiltoniar£3H. Hallberg et aLEJ studied the critical behavior of S" = 1/2 and S = 3/2 quantum spin chains 
^ I with periodic boundary conditions through extensive calculations of ground state correlation functions at different 
separations and different chain sizes L. Spin correlation functions in an open chain have also been calculated and 
compared_with results calculated from low-energy field theory, showing that estimates of the amplitudes can also be 
obtainedE^I. The approach describedUn this paper was applied to the spin-1/2 Heisenberg chain and a non-Hermitian 
supersymmctric (SUSY) spin chainE3. More recently, critical behavior of classical one-dimensional reaction-diffusion 
modelsO and the two-dimensional Potts modelH has been studied using the finite-size DMRG algorithm. Bulk 
and surface exponents of the Potts and Ising model have been obtained by using the DMRG to calculate correlation 
functions at different separations and collapsing curves obtained at different system sizest3. The SUSY chain describing 
the spin Quantum Hall effect (SQHE) plateau_transition was also examined in some detail. Critical exponents were 
extractedcj and compared to exact predictionscj. TherniOidja|iamic properties of other two-dimensional classical critical 
systems have also been studied by the DMRG methoc]c3~L3. Finally, Andersson et al. investigated the convergence 
of the DMRG in the thermodynamic limit for a gapless system of non-interacting fermionsc3. 

The method described in this paper combines the DMRG algorithm with finite-size scaling analysis, and yields 
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accurate critical exponents. The main advantage of the method is its simphcity. Only the calculation of ground state 
correlations near the rrviddle of chains with open boundary conditions are required. The relatively simple "infinite- 
size" DMRG algorithmB is particularly accurate for this job. In Sec. || we describe the method. The tight-binding 



model can be solved exactly and in Sec. Ill we use it to illustrate our scaling analysis. DMRG results are presented 
in Sec. ^ for the anisotropic S = 1/2 Heisenberg antiferromagnet and several critical exponents are obtained. An 
analytical calculation shows that multiplicative logarithmic corrections - which complicate the extraction of accurate 
critical exponents - may be avoided in some instances. In Sec. ^ the 5 = 1 antiferromagnetic spin chain is studied, 
focusing on the critical point that separates the Haldane and the dimerized phases. We conclude with a summary in 



Sec. VI, 



II. THE DMRG / FINITE-SIZE SCALING APPROACH 

We first describe how critical exponents may be obtained from a finite-size scaling analysis of chains with open or 
fixed boundary conditions. These boundary conditions are the simplest to implement in DMRG calculations. In the 
next subsection the DMRG algorithm itself is briefly described. 

A. Finite-Size Scaling 

To illustrate the sorts of power-law scaling we wish to examine, first consider the case of a spin chain with periodic 
boundary conditions that is at its critical point. The system can be moved away from criticality by turning on a 
uniform magnetic field, say in the x-direction, at each site: 

L 

HB = hY,S]. (1) 

This perturbation makes the correlation length finite: 

is cx \h\-^- . (2) 

Explicit dimerization, breaking the symmetry of translation by one site, also moves the system away from criticality. 
For a Heisenberg antiferromagnet, this can be realized by the addition of a staggering term R to the Hamiltonian: 

H=J2[l + {-iyR]SrS,+i. (3) 

The correlation length ^ in this case scales as 

c « \Rr • (4) 

Thus there are two independent exponents which correspond to these two perturbations of critical spin chains. Two- 
parameter scaling functions can be written for various observables and, for a finite system, these involve two dimen- 
sionless variables: the ratios L/^ and L/(^b- The induced dimerization, defined for now as the modulation of the x — x 
and y — y spin-spin correlations on even versus odd links, 

A = (-i)^- [{s^s^^, + sysy^^-{s^_,s^ + sy_,sy)] , (5) 

is of course independent of the site index for periodic chains, and scales as a function of the chain length L, the field 
/i, and the dimerization parameter R as: 

A{L, R, h) = sgn(i?) |i?r^ /A(i|i?r, L\h\-^) . (6) 

When the applied magnetic field is removed, h — Q, and this expression simplifies to: 

A(i,i?) = sgn(i?) |i?r- g^{L\Rn 

1"=^ R as i? -> 0, (7) 
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where the second hne follows from the fact that when the perturbation R is very small, or equivalently when the 
correlation length is larger than the system size, the net induced dimerization must be an analytic, linear, function 
of R. Therefore, for |a;| <C 1, the scaling function gA{x) is given by: 

9a{x) = Ixl-"-/"- ia,\x\'/'' + aM^'" + •■•); (8) 

the first term yields linear dependence of A in R in the R —* ^ limit, in agreement with Eq. ^ and the subsequent 
terms are higher order corrections. To recover the correct L-dependence, we must set 

= — - — ■ (9) 



The exponent xi\ and the correlation length exponent v satisfy the usual relation 

1 



(10) 



The applied magnetic field also polarizes the spins along the chain. The scaling form for the spin moment at each 
site is given by: 

(5-) =sgn(/i) \h\^- jB{L\RYMhD . (11) 
With no applied dimerization, i? = 0, and we expect the simple power-law: 

(5^) /i as /i . (12) 

Therefore, a;^ = (1 — aB)/vB- 

Alternatively, dimerization can be induced by open boundary conditions, and we take advantage of this fact to 
extract critical exponents. As depicted in Fig. open boundary conditions favor enhanced nearest-neighbor spin- 
spin correlations on the two outermost links. Chains of increasing length L = 4,6,8, .. . exhibit alternating patterns 
of dimerization on the interior bonds. Likewise, spin moments may be induced in the interior of the chain by applying 
a magnetic field to the ends of the chain. Strong applied edge magnetic fields completely polarize the end spins and 
induce non-zero and alternating spin moments along the chain. Alternatively, spin moments can be induced as before 
by a staggered magnetic field applied along the entire chain. Here however we consider only edge magnetic fields. 



• •^--^ • L = 4 

• • -• • L = 6 

• • -• • L = 8 
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FIG. 1. Extraction of critical behavior from finite-size effects. Dimerization of the nearest-neighbor spin-spin correlation 
function, indicated here by alternating strong (solid) and weak (dashed) bonds, is induced by the open boundary conditions. 
Spin moments are induced by the application of a magnetic field of strength ±/i to the two spins at the ends of the chains. 



We monitor the induced dimerization and spin moments at the center of the chain as the chain length L is enlarged 
via the DMRG algorithm. This scaling analysis is convenient because the relatively simple infinite-size DMRG 
algorithm applies to open chains and is most accurate at the center region of the chain where we focus our attentioe,. 
The induced dimerization and spin moments in the interior of the chain show power- law scaling at the critical point iS. 
Igloi and Rieger demonstrated power-law scaling for a variety of open boundary conditions (free, fixed and mixed)LZ. 
At the critical point i? = and h = Q the induced dimerization scales as a power-law with possible multiplicative 
logarithmic corrections: 

A(L/2) = L-"'^ (In L)y^ (a + A + . . . ) . (13) 

A similar expression holds for the induced spin moment at the center of the chain, {S^{L/2)), with the replacement 
of the exponents xa ^ xb and j/a ^ ys- 



B. Infinite-Size DMRG Algorithm 



The name "density-matrix renormalization-group" is something of a misnomer as the method is most accurate away 
from critical points, when there is an energy gap for excitations. It is helpful to think of the DMRG algorithm as a 
systematic variational approximation for the calculation of the ground state and/or low- lying excitations, principally 
in one dimension. The Hilbert space of a quantum chain generally grows exponentially with the chain length, and 
eventually must exceed available computer memory. The DMRG algorithm is an efficient way to truncate the Hilbert 
space; as the size of the space retained can be varied (up to machine limits) it is possible to ascertain the size of errors 
introduced by the truncation. _ 

For simplicity, we use the so-called "infinite-size" DMRG algorithmfi. As the algorithm has been described in some 
detail by White, we just sketch the essentials of the method. It begins with the (numerically exact) diagonalization 
of an open chain consisting of just four sites, each site having on-site Hilbert space of dimension D. For quantum 
spin chains D ^ 2S + 1, thus D = 2 for the spin-1/2 Heisenberg antiferromagnet. The chain is then cut through the 
middle into two pieces, one half of which is interpreted as the "system" and the other half as the "environment," the 
two parts combined being thought of as the entire "universe" of the problem, see Fig. ||. At this point the reduced 
density matrix for the system, of size DM x DM is constructed by performing a partial trace over the environment 
half of the chain. It is defined by: 

DM 

/'y (14) 

where ^^^z = are the real- valued matrix elements of the eigenstate of interest (the "target" which is often 

the ground state) projected onto a basis of states labeled by unprimed Roman index i which covers the system half 
of the chain and primed index i' which covers the environment half of the chain. The eigenvalues of the reduced 
density matrix are real, positive, and sum up to one; these are interpreted as probabilities. We keep only the M most 
probable eigenstates corresponding to the largest eigenvalues, and discard the remaining M{D — 1) eigenstates. The 
retained states form a new basis for the problem. Next, two new sites are added to the middle of the chain and the 
pieces are connected, yielding a chain of size L — 6. The process is then repeated by finding the targeted state of this 
chain, constructing the new reduced density matrix and again projecting onto the M most probable states. As the 
chain length grows in steps of two, the total Hilbert space dimension grows by a multiplicative factor of D^. None of 
the Hilbert space is thrown away until the chain grows large enough that its Hilbert space exceeds the space that is 
held in reserve, in other words until D^ > D^M^. The truncation process damages the outer regions of the chain the 
most, and the central region is treated most accurately. 
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FIG. 2. Schematic of one iteration in the infinite-size DMRG algorithm. At each step, a "new block" is formed for each half 
of the chain (the "old block" plus one additional site) and two more sites are added in the middle of the chain, increasing its 
length by two. 

One advantage of the method presented in this paper is that critical exponents are extracted from ground-state 
correlations only. Excited states are not needed for these exponents and there is no need to calculate the excitation 
gap. Furthermore, the finite size analysis described in the previous subsection takes advantage of the fact that the 
DMRG algorithm works best with open chains and treats the central region of the chain most accurately. The use of 
the more complicated finite-size algorithm might yield even more accurate results. However, we show below that we 
can calculate critical exponents to an accuracy of a few percent or better with the infinite-size algorithm. 



III. TIGHT BINDING MODEL AT HALF-FILLING 



As a simple first illustration of our finite-size scaling method we study the ordinary tight binding model of spinless 
fermions hopping from site-to-site along a chain at half-filling. Obviously, the DMRG algorithm is not needed in 
this case as we can solve the quadratic problem exactly via a Fourier transform. Due to particle-hole symmetry, at 
half- filling the chemical potential is zero. The correlation length exponent for this system is u = 1. A direct way to 
see this is by introducing the staggering parameter R to modulate the amplitude of the hopping matrix elements on 
even versus odd links: 

L-l 

H = t^[l + {-iyR]{c]cj+r+h.c.) . (15) 

3=0 

To diagonalizc the Hamiltonian, in the case of periodic boundary conditions cq = Cl, we introduce separate fermion 
operators for even and odd sites as follows: 

C2J-1 = e2j (16) 
After the Fourier transformation to momentum-space, the Hamiltonian can be written as: 
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i? = t^|[(l-i?) + e'*'(l + ^)] diek + [{l-R) + e-^'''{l + R)] e^dfej, (17) 
where the lattice spacing a — 1. For each k, diagonaUzation of the 2x2 matrix yields the dispersion relation: 



Cfe = ±2i^l - (1 - i?2)sin^(fc). (18) 

At half-filling the ground state has all states with < occupied. The left and right Fermi points are, respectively, 
kp = ±7r/2. Hence the gap m = 2t|i?|. As the correlation length ^ c>c oc \R\~^ we obtain = 1. Since 

ly^^ = 2 — a; A — 1, the dimerization exponent xa = 1- 

We now reproduce this result using the finite-size scaling method applied to open chains. We consider a finite chain 
of length L with open boundary conditions and calculate the induced dimerization A(j) = (—I)-' {cjCj+i — 
around the chain center j ~ L/2, and extract its leading dependence on L. Open boundary conditions are imposed 
by using the Fourier transform 



1 



L 



=1 



km = -i—rm, m = l,2,...,L (19) 

as this enforces co = cl+i = 0. Filling all of the negative energy states at half-filling, the expectation value of the 
dimerization at L/2 can be found by straightforward calculation: 

1 ^ 

A(L/2) oc j-^ [ cos[fc™(L + 3)] - cos[A:,„(i + 1)]] . (20) 

This sum can be evaluated numerically with the result that ^ 1 as L ^ 00 as shown in Fig. ^, in agreement with 
the explicit calculation for the periodic chain. It is also easy to show that open chains with an odd number of sites 
have vanishing induced dimerization at the center of the chain, as expected by the symmetry of refiection about the 
central site. 




FIG. 3. Induced dimerization in the tight binding model. A log-log plot of the dimerization at the center of the chain (Eq. 
po| ) is shown as a function of the chain length L. For more complicated systems, the DMRG algorithm is employed to calculate 
A(L/2) numerically. In the inset, the dimerization exponent calculated from the slope of the curve shown in the main graph 
is plotted as function of l/ln(L). For small L there are subleading corrections to scaling, but xa — > 1.0 as the chain length 
increases. 
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The induced density moment can likewise be obtained either directly by studying the effects of a staggered chemical 
potential /istag (which doubles the size of the unit cell from one to two sites and thus generates a gap m = 
or by the inclusion a local chemical potential fj, at the two ends of the chain: 

H^H-fi (cjco + c{_,CL-i) . (21) 

Again the system consists of L sites, the site index running from to L — 1, and there are open boundary condition 
at J = and j = L — 1. For large /i ^ 0, the boundary condition is equivalent to enforcing unit occupancy at the 
chain ends, no = ^l-i = 1- This boundary condition is satisfied by the Fourier transform 



^ / ik„ 



(22) 



with 



(23) 



Again it is a simple exercise to calculate the occupancies. At the chain ends we obtain: (cqCq) = {c\^iC'L-i) = 1 in 
agreement with the boundary condition. At the center of the chain the occupancy can be evaluated analytically, 

L-1 

(ct(i/2)c(L/2)) = j—^ ^ [1 + cos(fc™L)] . (24) 

It scales as {c'{L/2)c{L/2)) — 1/2 c>c L^^. Hence = a^s = 1 in agreement with the direct calculation of these 
exponents. 



IV. SPIN-1/2 ANTIFERROMAGNET 



We next turn to the study of a richer system: spin-1/2 antiferromagnetic chains. We begin with the XY model, 
which can be solved exactly by a Jordan- Wigner mapping to the tight binding model. We then study the anisotropic 
XXZ model. The isotropic Heisenberg model is treated separately as there are complicating multiplicative logarithmic 
corrections to scaling at the isotropic point. 



A. XY model 



The Hamiltonian for the spin-1/2 XY model. 



L-2 



qV qV 1 



(25) 



can be writte n in terms of spinless fermion creation and annihilation operators and cj via the Jordan- Wigner 
transformationcS. An up spin in the z-direction at site i then corresponds to having the site occupied by a fermion, 
while spin down corresponds to an empty site. The Hamiltonian of Eq. |2^ is mapped to a nearest-neighbor tight 
binding Hamiltonian with t = J. Based on our analysis in the previous section we can conclude that v — I for the 
XY model. 

Fig. ^ presents our DMRG results for the induced dimerization and induced spin moments, in the x- and in the 
z-directions, at the center of the chain as a function of the chain length, L. 
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FIG. 4. Log-log plots of induced dimerization and induced spin moments in the z-direction and x-direction (with the edge 
magnetic field applied in the z- and in the x-directions, respectively) at the center of the chain for the spin-1/2 XY model. The 
magnitude of the applied edge magnetic field is h = 1.0 in both cases and the number of block states kept is M — 128. 

The exponents are obtained from the slopes of the curves shown in Fig. ^. The induced dimerization exponent 
for A(L/2) is close to 1 {xa = 0.99 ± 0.01) as expected from the relation — 1/(2 — xa)- The slope of the log-log 
plot of the induced spin moment in the z-direction is also close to 1 (xs = 1.01 ± 0.02). This result is also expected 
since it is equivalent to the exponent for the induced density moment in the tight binding model as discussed in the 
previous section. In the case of the induced spin moment in the x-direction, the exponent is 0.248 ± 0.003. This value 
compares well with the exact number of 1 /4 as derived in the next section. 



B. XXZ model 



Next consider the nearest-neighbor, spin-1/2 XXZ Heisenberg antiferromagnet: 



H = J 



L-2 

E 



[S^S: 



3 



qy qy 



i+iJ 



(26) 



Anisotropy in the coupling between the z-components of the spins may be varied by chan; 
Jordan- Wigner transformation, the XY terms again yield the tight binding Hamiltonian. 
therefore occur near the two Fermi points at A: = ±7r/2a. We may treat the non-Gaussian 7 
and focus on excitations around these Fermi points by defining left and right moving low-energ; 
the continuum limit and keeping only the low-energy modes, the tight binding term is then 
the massless fermions. The S^S^^i term is quartic in the fermion operators. Integrating out 
it will renormalize the fermion velocity and also contain interaction terms. 

We then implement Abelian bosonization, with UV cutoff a. The effective Hamiltonian is a sine-Gordon model (a 
derivation can be found in Ref. 49): 



gmg 7. Performing the 
Low-energy excitations 
term as a perturbation 
J quasiparticles. Taking 
effectively described by 
the high-energy modes. 



H = Hn- 



27ra2 



dx cos\v87r(h(x) 



where 



Hq = u I dx 



[OA? 

K 



(27) 



(28) 
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Here u = 2Ja = 2a is the bare Fermi velocity and the constant K = 1 + yo/2 depends on the anisotropy 7. The XY 
hmit correspond to = 0. 

The long distance behavior of the staggered part of and are given in terms of the boson fields as: 



S'^ix) « {~lf/°'cos[(j)ix)/R] 

S~{x)^{-lf^°'e'^''^'^^'='> (29) 



where the radius R is given b; 



First consider the anisotropic case 7 7^ 1. The isotropic case has logarithmic corrections to scaling that are dealt with 
in the next section. For 7 > 1 the interaction term is relevant and the system is gapped, and in the Ising universality 
class. Indeed, in the limit 7 — > cx) it is the Ising model. For 7 < 1 the interaction term is irrelevant, the system is 
gapless and A(L/2) and {S^{L/2)) should exhibit power law decay, with no log corrections as there are no marginal 
operators. The log- log plots of Fig. |5| (a) show the induced spin moment in the x-direction at the chain center for 
different values of the anisotropy 7. The edge magnetic field in the x-direction is fixed, h — 1.0. As expected, for 
7 > 1 there is exponential decay and in the cases 7 < 1 the exponents xb{'^) are found by fitting the curves in Fig. 
^ (a) to the form of Eq. |l^. The exponents ys are set equal to zero, the higher order corrections are included and 
give very small deviations from a simple Uiiear fit. In Fig. ^ (b) the exponents xsij) are compared to the exact 
value xb{'~i) = itR'^{'j) obtained by AffleckEj. Agreement is found at the percent level. Affleck derived the exponent 
as follows. The edge magnetic field in the x-direction applied at j = corresponds to a term 

Hb = -hS^'iO) = -constant x h cos[V2^4>{0)] (31) 

in the Hamiltonian. For sufficiently large h the energy is minimized by setting 

^(O) = O^0fl(O) = (/)L(O) . (32) 

Regarding as an analytic continuation oi (pL, we may identify 

<I>b{x) ^ <f>L{-x) . (33) 

Using this boundary condition, the induced spin moment is given by 

(5-(j)) « (_i)./"(e^2.fl0.o-)^-.2.fl0.(-,)^ ^ ^^yTi,) ■ (34) 
For the XY model (7 = 0), the induced spin moment in the x-direction therefore decays with exponent 7ri?'^(0) — 1/4. 
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FIG. 5. Spin-1/2 XXZ Heisenberg antiferromagnetic chain, (a) DMRG results (Af = 128 and h = 1.0) for the induced 
spin moment in the x-direction at center of chain for different values of anisotropy 7, showing power law decay for 7 < 1 and 
exponential decay for 7 > 1 and (b) corresponding exponent 2:3(7) for 7 < 1- The solid line is the exact value -kR (7). The 
exponents are obtained by fitting the curves in (a) to the form of Eq. ^ with yB = since there are no logarithmic corrections 
away from the isotropic point. 

A log-log plot of the induced dimerization at the center of the chain for various values of the anisotropy 7 is shown 
in Fig. ^ (a). The free boundary condition at the chain ends corresponds to setting: 

So = Sl+1 = . (35) 
This condition translates to 4>r{x) — —(^^{—x) + nR in terms of the boson fields0 which yields: 
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In Fig. H (b) the exponents obtained from the slopes of the curves in Fig. ^ (a) are plotted against the exact values 
XA = l/47ri?^(7). Again agreement is found at the percent level. 
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FIG. 6. Spin-1/2 XXZ Heisenberg antiferromagnetic chain, (a) DMRG results (M — 128) for the induced dimerization at 
center of chain for different values of anisotropy 7 and (b) corresponding exponent a;A(7), obtained by fitting to Eq. ^ with 
yA = 0. The solid line is the exact value 1/4:-kR'^{j). 



Another quantity of interest is the sum, instead of the difference, of the spin-spin correlation function on adjacent 
bonds near the center of the chain: 
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e{L/2) ^ - ( {Sl/2 ■ Sl/2+i) + {Sl/2^iSl/2) ) , 



(37) 



which at 7 = 1 equals the energy density per bond and therefore does not vanish in the thermodynamic hmit. Fig. |^ 
is a plot of e(L/2) as a function of the system size L at the isotropic point 7=1. 
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FIG. 7. Sum of the two central bonds for the spin-1/2 XXZ Heisenberg antiferromagnetic chain at the isotropic point 7 : 
The block size is M = 256. 



As expected, this quantity approaches a constant value e(oo) in the thermodynamic limit. After subtracting the 
extrapolated value at L ^ 00, e.{L/2) too exhibits power law decay of the form of Eq. |3|. The constant e(oo) can 
be found by an iteration process. Starting with an initial value for e(oo) obtained from a rough extrapolation of the 
curve in Fig. 0, we fit the subtracted value e(L/2) — e(oo) to a power-law form. The extrapolated value e(oo) is 
then adjusted slightly until an optimal fit to a pure power law is attained. The extrapolated value found this way is 
e(oo) = -0.443148 and Fig. | shows the power law behavior of the subtracted quantity. 
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FIG. 8. Log-log plot of e(L/2) after subtracting the extrapolated value e{oo) = —0.443148. The fit is to a straight line of 
slope 2.14. 

We obtain an exponent of 2.1 ± 0.1 in the scaling of e{L/2) — e(oo). This is as expected from the linear dispersion 
relation of Heisenberg antiferromagnets: in a Lorentz-invariant theory the energy density operator has dimension 2. 

The DMRG result for the eneipgy per bond is extremely accurate and can be compared with the exact value obtained 
from the Bethe ansatz solutionEj of e = 1/4 — ln2 = —0.44314718. It is crucial to note that the open boundary 
conditions induce staggering in the strength of the bonds along the chain. To eliminate this effect, the energy per 
bond must be calculated as the average of the bond energy from two consecutive bonds at the cepter of the chain. 
Suggestions that infinite-size DMRG results for the center region of the chain are not very accurately appear to have 
failed to take this effect into account. We have also checked our results at different anisotropies. For the XY case 
(7 = 0.0), we obtain e(oo) = —0.318310 extrapolating from chains up to L = 200 and M = 128 and the exact resultEHl 
is -I/tt = -0.3183099. 



C. Logarithmic Corrections to Scaling 

In the isotropic XXX limit, the interaction cos[\/8TT(f>{x)] in the low-energy effective Hamiltonian Eq. |2^ becomes 
marginal and can generate multiplicative logarithmic corrections to scaling. In this section we calculate its effect on 
the scaling of the induced spin moment {S^{L/2)) when an edge magnetic field Hb in the x-direction is applied. Can- 
cellations occur and in this case there are no multiplicative ln(L) corrections. As a practical matter, the cancellation 
of the logarithmic corrections means that numerical calculations of the exponent xb are particularly precise. We note 
that finite-size scaling of thepSpia-spin correlation function has been previously calculated for a spin- 1/2 chain with 
periodic boundary conditionaEJo. 

The coupling constants in the sine-Gordon Hamiltonian (Eq.j-g^) rcnormalizc under a change of the ultraviolet 
cutoff a — > ae' according to the renormalization group cquationsES"^ 

^ = -yS)v^{i)- (38) 

As noted in the previous section, a large edge magnetic field applied at a: = in the x-direction enforces the boundary 



condition (j)R{x) — (j)Li~x) (Eq. |33|). Thus 



(^"(x)) - (-l)"/''(cos[\/2^0(a;)]) - (^e'^'^^^^^e-'^'^^^-^^). (39) 
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For the free theory, which corresponds to the XY model = 0, the induced spin moment is simply 

{S^x))o^exp[-KUL{2x)] (40) 

where 

C/^(x) = iln(^). (41) 
2 a 

But in the general XXZ case we ascectain the effect of the marginal operator by following a procedure similar to one 
developed by Giamarchi and SchulzE3 who calculated correlation functions for finite periodic chains. We first define 
the function: 

F{x) = e^^^'^^^\S^ix)). (42) 

At the XY point = Q clearly F[x) — 1. For small x, an expansion of F in powers of converges, and for sufficiently 
small coupling y^, F{x) ^ 1. Upon rescaling, the function F{x) also depends on the new length scala_and on the 
rescaled coupling constants yo(0 E^nd y,f,{l). By an argument similar to the one employed by KosterlitzEll, the effect 
of rescaling a ^ e' a is: 

F{x,ae',y(l)) = I(dl,y{l))F{x,ae'+''',y{l + dl)) , (43) 

where y{l) denotes all the couplings as function of the scaling parameter I. The rescaled short distance cutoff is then 
a{l) = e' a, where a is the initial cutoff. Rescaling can be repeated until a{l) ~ x, at which point we have: 

F{x,x,yiHx/a))) = Oil) . (44) 

The contributions to the function F from repeated rescalings, until a(^) reaches x, can be written explicitly as: 

l=\n{x/a) , Ani'x/a 



Fix,a,y{a))= Hdl,y{l)) ^ exp \ [ In [/(d/, y(0)] < 

1=0 ^ 



dl )■ . (45) 



We proceed to calculate the function /. First we expand {S^{x)) in powers of y^, writing it in terms of averages 
with respect to the free Hamiltonian, 

(5"(x)) ~ e-^^^(2.) + / d^x'{S'={x)cos[V8^c^{x')])o + 

1/ y4> \2 

The averages are given by 

(S"^(x)cos[\/8^<?!)(a;')])o = 

and 



d^xi / d^X2{S'={x) cos[V87r0(a;i)] cos[V87r(/)(a;2)])o + ■ • ■ (46) 



(S'^ix) cos[V87r(/)(xi)]cos[V87r0(a;2)])o 

i exp [ - KUl{2x) + AKUl{2xi) + AKUl{2x2) - AKU{xi + X2) - AKU{xi ~ X2)] (47) 

The 0{y^) term can be simplified by assuming that the main contribution conies from configurations where xi and 
X2 are very close to each otheril'0. Introducing new integration variables 

r = xi — X2 

R^'-^ (48) 
and expanding C/l in powers of r, which is assumed to be small, 

Ul{2xi) = ULi2R + r) = Ul{2R) + r Or Ul{2R) + . . . , 

Ul{2x2) = UL{2R-r) = Ul{2R) - r Or Ul{2R) + . . . , (49) 
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we obtain the average 



{S^'ix) cos[V8^<j){xi)] cos[\/8^(/i(x2)])o ^ exp [~KUl{2x) - 4:KU{r)] 



The dependence on R cancels out. The expansion Eq. 



(S^ix)) 



I becomes 



(50) 



(51) 



where = / dR is a measure of the linear size of the system. Next consider the effect of rescaling a' — ae'^', where 
dl is infinitesimal. Using 



dx 



dx 



dx 



we obtain: 



4a' 



Matching this result with Eq. |43|, we find: 

I{dl,yo{l),y^{l)) ~ exp 



4a2 



dl 



hence from Eq. 45 and Eq. k2, we have 



(5"(x)) - exp <{ - KUl{2x) + J -^dl 



4a2 



Using the RG equations (Eq. p^), the solution at large I is y^{l) ^ l/l and 



{S^{x)) 



'1/2 



exp 



ln{x/a} -, 

d/[0(-) 



(52) 



(53) 



(54) 



(55) 



(56) 



There are no multiplicative ln(a;) corrections, as these would require terms of order 0{l/l) in the integrand inside the 
exponential in Eq. [s^. In our calculation, 0{l/l) terms do not appear, only 0{1/P) and higher-order terms. As a 
check, we can repeat the same procedure for {S^{x)), with the edge field now oriented in the z-direction. Of course 
this should give the same result since the system is isotropic, but as the Jordan- Wigner transformation picks the 
z-direction as the spin quantization axis, the equivalence is not obvious, and the check is non-trivial. Again, explicit 
calculation shows that 0(1//) terms do not arise. This result is in reasonable agreement with our numerical results. 
Fitting the DMRG data (see Fig. |) to the form Eq. |l|, we obtain xb — 0.485 ± 0.01 with a small non-zero value for 
the log exponent ys = 0.06 ± 0.01. By contrast, in the case of the induced dimerization we obtain xa = 0.57 ± 0.01 
and t/A = 0.10 ± 0.05. The error was estimated from deviations obtained by fitting the M — 128 data over different 
ranges of L (4 < i < 600) and by comparison with M = 64 data (4 < L < 300). Results are systematically improved 
by increasing the value of M. 
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FIG. 9. Log-log plot of the induced spin moment (h = 1.0) and induced dimerization at center of chain for the isotropic 
S — 1/2 Heisenberg antiferromagnet. Data obtained using block Hilbert space size of Ai — 128 and M = 64 are plotted for 
comparison. Note the slight curvature in A(L/2) which signals the existence of multiplicative logarithmic corrections to scaling. 
The log is absent in the plot of {S^{L/2)). The discrepancy in the induced dimerization at large chain length is due to the 
truncation of the block Hilbert space. 

Finite-size scaling behavior for the XXX model with open boundary condition^ and periodic boundary conditions0 
were obtained from DMRG calculations of ground state energies and correlation functions {S^{x)S^{x+r)) for different 
system sizes and separations r. In our approach, critical exponents are extracted from expectation values at the center 
of the chain only. The chain size is increased via the infinite-size DMRG method. It is also advantageous to extract 
power law exponents when there are no logarithmic corrections. 



D. Conformal Anomaly 

Finally, we may calculate the value of the conformal anomaly, c. We note that the central charge of the RSOS 
modeCj and of the spin-3/2 Heisenberg chaino, which is in the same universality class as the spin-1/2 chain, have 
previously been obtained using the DMRG. The conformal anomaly can be extracted by finding the coefficient of the 
1/L finite-size correction to the free energy, equivalent at zero temperature to the ground state energy. We fit the 
ground state energy Eq{L) to the following form: 

Eq{L)^AL + B + C/L + ... (57) 

The extensive contribution, proportional to A, and the constant term B are non-universal. At the isotropic point 7 = 1 
our results for the case of blocksize M — 128 and for chain lengths in the range 30 < L < 100 yield C « —0.323. To 
relate this coefficient to the conformal anomaly wc must normalize it by dividing by the speed of low- lying excitations, 
V. The speed can be obtained by extrapolation to the thermodynamic limit of the gap to the lowest-lying excitation 
multiplied by the chain length: 

. = hm ^^m^^ . (58) 

L — ^00 TT 

We find v — 2.44. Now for open boundary conditions^, 

-24C 



1.01 . (59) 
nv 

This compares well with the value of c = 1 appropriate for a single boson or the pair of left and right moving fermions. 



16 



V. SPIN-1 CHAIN 



As a final example we apply the DMRG / finite-size scaling method to the isotropic spin-1 antiferromagnetic chain. 
This problem is more challenging numerically as the on-site Hilbert space now has dimension D = 3 instead oi D — 2. 
The most general nearest- neighbor Hamiltonian for the spin-1 chain includes the possibility of a biquadratic spin-spin 
interaction term: 

L-2 

H =^[cos9 Sj ■ Sj+i + sin {Sj ■ Sj+i)^] (60) 
j=o 

The phase diagram can be represented on a circle parameterized by 6 as depicted in Fig. |l^. 




FIG. 10. Phase diagram for isotropic, nearest-neighbor, S=l chain, parameterized by the angle 9. A critical point at 
9 = — 7r/4 separates the Haldane and dimerized phases, both of which are gapped. At 9 = tan~^(l/3), the ground state is a 
valence bond solid (VBS), as depicted in the schematic. Each oval represents an atom; the two black dots inside the oval are 
two electrons symmetrized into a triplet 5 = 1 spin state. Rectangles represent singlet bonds which encompass two electrons 
on adjacent sites. In the dimerized phase, singlet correlations are instead enhanced on alternating dimers. 

Generically there-is a gap to excitations in the antiferromagnetic region of the phase diagram, in accord with the 
Haldane conjectureEa. The point 9 = corresponds to the usual pure bilinear Heisenberg antiferromagnet. At the 
point tanO = 1/3 the Hamiltonian can be written a sum of positive-definite projection operators, and the exact 
ground state is the AKLT valence bond solid (VBS)!^. Negative sinO favors dimerization, as the energy is minimized 
by concentrating singlet correlations on isolated dimers. The dimerized phase also is gapped: a dimer must be broken 
to generate a spin excitation. The point, that separates the dimerized and Haldane phases lies at 9 = — 7r/4 and can 
be solved exactly by the Bethe ansatzL^Ej. The chain is quantum critical at this integrable point. The ground state 
is non-degenerate here as well as in the dimerized and Haldane phases. 

DMRG calculations clearly delineate the two massive phases and the critical point separating them, even for 
relatively small block Hilbert sizes M. In Figs. |l^ and ^ the blocksize M = 81 for the massive phases. Thus the 
results are numerically exact up only to chain lengths L = 10. For chain lengths L > 10 the Hilbert space is truncated 
via the DMRG algorithm. To increase accuracy, results at the critical point were obtained with a larger Hilbert size 
for the blocks, M = 256. 



17 



1.000 r 




L 

FIG. 11. Induced spin moment at center of a spin-1 chain for different values of 6. Here ft = 1. Note the power law decay at 
the critical point (9 = — 7r/4). Exponential decay occurs at the AKLT valence bond solid point [6 = arctan(l/3) = 0.321751] 
within the Haldane phase and in the dimerized phase at = — 7r/2. Results at the critical point were obtained with M = 256, 
while M = 81 for the other two cases. 




L 

FIG. 12. Induced dimerization at center of a spin-1 chain. Again we set M = 256 at the critical point and M = 81 for the 
two massive phases. There is power law decay at the critical point {9 = — 7r/4) and exponential decay at the AKLT point. As 
expected the dimerization approaches a non-zero constant in the dimerized phase (9 = — 7r/2). 

The induced spin moment at the center of the chain decays exponentially in both the Haldane and the dimerized 
phases, as expected. The induced dimerization at the chain center also decays exponentially in the Haldane phase. 
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but approaches a non-zero constant in the dimerized phase as it must. Power law decay in both observables occurs at 
the critical point. Fitting the M = 256 data shown in Fig. [ij at the critical point 9 = — 7r/4 we obtain dimerization 
exponents — 0.37 ± 0.01 and j/a = 0.3 ± 0.05, reflecting the apparent presence of a marginal interaction and 
consequent multiplicative logarithmic corrections to scaling. Likewise, for a field of ft, = 1.0 applied to the chain ends, 
the exponents for thj-spin operator are xb — 0.34 ±0.01 and ys — 0.23 ±0.05. The values of the exponents compare 
to the exact valueg£jO xa = 3/8 ~ 0.375 and xb = 3/8. To the best of our knowledge there are no analytic results 
at the integrable point 9 = — 7r/4 (which corresponds to a fc = 2 SU(2) WZW model) on the size of the logarithmic 
corrections j/a and ys, at least for open boundary conditions. 

Finally, we may repeat the analysis of the conformal anomaly described above in subsection IV D| for the case of 
the spin-1 chain at its critical point. For M = 256 and fitting over chain lengths 10 < L < 26 we find that the speed 
of excitations is w = 3.69, C = 0.508, and hence c — 1.05. This value is close to its exact value of 1, demonstrating 
that the conformal anomaly can be reliably extracted even from relatively short chains. 



VI. CONCLUSION 



We have presented a simple method for studying critical behavior of quantum spin chains. Accurate critical 
exponents can be extracted. For small on-site Hilbert space sizes {D — 2 for the spin- 1/2 chain and _D = 3 for 
spin-1 chains) the method does not require supercomputers. Results can be systematically improved by increasing 
the size of M , the dimension of the Hilbert space retained in the blocks, up to limits set by machine memory and 
speed. The DMRG method works best for massive, non-critical, systems, but it is also quite accurate even at critical 
points. Critical exponents can be calculated at percent level accuracy. We showed that the leading multiplicative 
logarithmic correction to the scaling of the induced spin moment cancels out in the case of the isotropic spin-1/2 
Heisenberg antiferromagnet. Thus accurate exponents can sometimes be found numerically despite the presence of 
marginal interactions. 

Use of the "finite-size" DMRG algorithm might improve the method, but good results were obtained with the 
relatively simpler "infinite-size" DMRG algorithm. The reason for this is that the finite-size scaling method employed 
here focuses on the scaling of observables near the center of the chain only, where the "infinite-size" algorithm is 
particularly accurate. The method can be used to study new systems. For example, several non-interacting but 
disordered electjon systems, like the integer and spin quantum Hall transitions, can be described by supersymmetric 
HamiltoniansEllH. In a paper which foUowsa, we employ the combined DMRG/finite-size method in combination 
with analytic calculations to understand the behavior of these supersymmetric spin chains. 
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